The Globular Cluster System in the Inner Region of the Giant Elliptical Galaxy NGC 4472

We present a study of globular clusters in the inner region of the giant elliptical galaxy NGC 4472, based on the HST WFPC2 archive data. We have found about 1560 globular cluster candidates at the galactocentric radius r < 4 arcmin. V-(V-I) diagram of these objects shows a dominant vertical structure which consists obviously of two components: blue globular clusters (BGCs) and red globular clusters (RGCs). The luminosity function of the globular clusters is derived to have a peak at V(max)=23.50+/-0.16 from Gaussian fitting. The distance to NGC 4472 is estimated to be d=14.7+/-1.3 Mpc.The peak luminosity for the RGCs is similar to that for the BGCs, which indicates that the RGCs may be several Gyrs younger than the BGCs. The mean luminosity of the bright BGCs decreases by 0.2 mag with increasing galactocentric radius over the range of 9 arcmin, while that of the RGCs does not. The observed color distribution of these globular clusters is distinctively bimodal with peaks at (V-I) = 0.98 and 1.23. The mean observed color of all the globular clusters with V < 23.9 mag is derived to be (V-I)=1.11. These colors are exactly the same as those for the globular clusters in M87. It is found that the relative number of the BGCs to the RGCs is increasing with the increasing galactocentric radius. Surface number density profiles of both the BGCs and RGCs get flat in the central region, and the core radii of the globular cluster systems are measured to be r_c = 1.9 arcmin for the BGCs, r_c = 1.2 arcmin for the RGCs, and r_c = 1.3 arcmin for the total sample, which are much larger than the stellar core of the galaxy. In general the properties of the globular clusters in the inner region of NGC 4472 are consistent with those of the globular clusters in the outer region of NGC 4472.Comment: 27 pages (AASLaTex), 22 Postscript Figures, Accepted for Publication in the Astronomical Journal, Jul. 31st, 200

On the Self-Consistent Response of Stellar Systems to Gravitational Shocks

We study the reaction of a globular star cluster to a time-varying tidal perturbation (gravitational shock) using self-consistent N-body simulations and address two questions. First, to what extent is the cluster interior protected by adiabatic invariants. Second, how much further energy change does the postshock evolution of the cluster potential produce and how much does it affect the dispersion of stellar energies. We introduce the adiabatic correction} as ratio of the energy change, , to its value in the impulse approximation. When the potential is kept fixed, the numerical results for the adiabatic correction for stars with orbital frequency \omega can be approximated as (1 + \omega^2 \tau^2)^{-\gamma}. For shocks with the characteristic duration of the order the half-mass dynamical time of the cluster, \tau < t_{dyn,h}, the exponent \gamma = 5/2. For more prolonged shocks, \tau > 4 t_{dyn,h}, the adiabatic correction is shallower, \gamma = 3/2. When we allow for self-gravity and potential oscillations which follow the shock, the energy of stars in the core changes significantly, while the total energy of the system is conserved. Paradoxically, the postshock potential fluctuations reduce the total amount of energy dispersion, . The effect is small but real and is due to the postshock energy change being statistically anti-correlated with the shock induced heating. These results are to be applied to Fokker-Planck models of the evolution of globular clusters.Comment: 20 pages; ApJ 513 (in press

The Globular Cluster System of the Spiral Galaxy NGC7814

We present the results of a wide-field photometric study of the globular cluster (GC) system of the edge-on Sab spiral NGC7814. This is the first spiral to be fully analyzed from our survey of the GC systems of a large sample of galaxies beyond the Local Group. NGC7814 is of particular interest because a previous study estimated that it has 500-1000 GCs, giving it the largest specific frequency (S_N) known for a spiral. Understanding this galaxy's GC system is important in terms of our understanding of the GC populations of spirals in general and has implications for the formation of massive galaxies. We observed the galaxy in BVR filters with the WIYN 3.5-m telescope, and used image classification and three-color photometry to select GC candidates. We also analyzed archival HST WFPC2 images of NGC7814, both to help quantify the contamination level of the WIYN GC candidate list and to detect GCs in the inner part of the galaxy halo. Combining HST data with high-quality ground-based images allows us to trace the entire radial extent of this galaxy's GC system and determine the total number of GCs directly through observation. We find that rather than being an especially high-S_N spiral, NGC7814 has <200 GCs and S_N ~ 1, making it comparable to the two most well-studied spirals, the Milky Way and M31. We explore the implications of these results for models of the formation of galaxies and their GC systems. The initial results from our survey suggest that the GC systems of typical ellipticals can be accounted for by the merger of two or more spirals, but that for highly-luminous ellipticals, additional physical processes may be needed.Comment: 28 pages, incl. 4 figures; accepted for publication in The Astronomical Journal, November 2003 issu

Dynamical Evolution of Globular Clusters in Hierarchical Cosmology

We probe the evolution of globular clusters that could form in giant molecular clouds within high-redshift galaxies. Numerical simulations demonstrate that the large and dense enough gas clouds assemble naturally in current hierarchical models of galaxy formation. These clouds are enriched with heavy elements from earlier stars and could produce star clusters in a similar way to nearby molecular clouds. The masses and sizes of the model clusters are in excellent agreement with the observations of young massive clusters. Do these model clusters evolve into globular clusters that we see in our and external galaxies? In order to study their dynamical evolution, we calculate the orbits of model clusters using the outputs of the cosmological simulation of a Milky Way-sized galaxy. We find that at present the orbits are isotropic in the inner 50 kpc of the Galaxy and preferentially radial at larger distances. All clusters located outside 10 kpc from the center formed in the now-disrupted satellite galaxies. The spatial distribution of model clusters is spheroidal, with a power-law density profile consistent with observations. The combination of two-body scattering, tidal shocks, and stellar evolution results in the evolution of the cluster mass function from an initial power law to the observed log-normal distribution.Comment: 5 pages, proceedings of IAU 246 "Dynamical Evolution of Dense Stellar Systems", eds. Vesperini, Giersz, Sill

Formation and evolution of clumpy tidal tails around globular clusters

We present some results of numerical simulations of a globular cluster orbiting in the central region of a triaxial galaxy on a set of 'loop' orbits. Tails start forming after about a quarter of the globular cluster orbital period and develop, in most cases, along the cluster orbit, showing clumpy substructures as observed, for example, in Palomar 5. If completely detectable, clumps can contain about 7,000 solar masses each, i.e. about 10% of the cluster mass at that epoch. The morphology of tails and clumps and the kinematical properties of stars in the tails are studied and compared with available observational data. Our finding is that the stellar velocity dispersion tends to level off at large radii, in agreement to that found for M15 and Omega Centauri.Comment: LaTeX 2e, uses AASTeX v5.x, 40 pages with 18 figures. Submitted to The Astronomical Journa

The growth of galaxies in cosmological simulations of structure formation

We use hydrodynamic simulations to examine how the baryonic components of galaxies are assembled, focusing on the relative importance of mergers and smooth accretion in the formation of ~L_* systems. In our primary simulation, which models a (50\hmpc)^3 comoving volume of a Lambda-dominated cold dark matter universe, the space density of objects at our (64-particle) baryon mass resolution threshold, M_c=5.4e10 M_sun, corresponds to that of observed galaxies with L~L_*/4. Galaxies above this threshold gain most of their mass by accretion rather than by mergers. At the redshift of peak mass growth, z~2, accretion dominates over merging by about 4:1. The mean accretion rate per galaxy declines from ~40 M_sun/yr at z=2 to ~10 M_sun/yr at z=0, while the merging rate peaks later (z~1) and declines more slowly, so by z=0 the ratio is about 2:1. We cannot distinguish truly smooth accretion from merging with objects below our mass resolution threshold, but extrapolating our measured mass spectrum of merging objects, dP/dM ~ M^a with a ~ -1, implies that sub-resolution mergers would add relatively little mass. The global star formation history in these simulations tracks the mass accretion rate rather than the merger rate. At low redshift, destruction of galaxies by mergers is approximately balanced by the growth of new systems, so the comoving space density of resolved galaxies stays nearly constant despite significant mass evolution at the galaxy-by-galaxy level. The predicted merger rate at z<~1 agrees with recent estimates from close pairs in the CFRS and CNOC2 redshift surveys.Comment: Submitted to ApJ, 35 pp including 15 fig

Determining the galactic mass distribution using tidal streams from globular clusters

We discuss how to use tidal streams from globular clusters to measure the mass distribution of the Milky Way. Recent proper motion determinations for globular clusters from plate measurements and Hipparcos astrometry provide several good candidates for Galactic mass determinations in the intermediate halo, far above the Galactic disk, including Pal 5, NGC 4147, NGC 5024 (M53) and NGC 5466; the remaining Hipparcos clusters provide candidates for measurements several kpc above and below the disk. These clusters will help determine the profile and shape of the inner halo. To aid this effort, we present two methods of mass determination: one, a generalization of rotation-curve mass measurements, which gives the mass and potential from complete position-velocity observations for stream stars; and another using a simple chi^2 estimator, which can be used when only projected positions and radial velocities are known for stream stars. We illustrate the use of the latter method using simulated tidal streams from Pal 5 and find that fairly accurate mass determinations are possible even for relatively poor data sets. Follow-up observations of clusters with proper motion determinations may reveal tidal streams; obtaining radial velocity measurements would enable accurate measurements of the mass distribution in the inner Galaxy.Comment: 21 pages, 6 figures, published in A